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ABSTRACT 



We present Keck II spectroscopy of an optical mHz quasi-periodic oscillations 
(QPOs) in the lightcurve of the X-ray pulsar binary Hercules X-l. In the power spec- 
trum it appears as 'peaked noise', with a coherency ^2, a central frequency of 35 mHz 
and a peak-to-peak amplitude of 5%. However, the dynamic power spectrum shows it 
to be an intermittent QPO, with a lifetime of ^hundred seconds, as expected if the 
lifetime of the orbiting material is equal to the thermal timescale of the inner disk. 
We have decomposed the spectral time series into constant and variable components 
and used blackbody fits to the resulting spectra to characterise the spectrum of the 
QPO variability and constrain possible production sites. We find that the spectrum 
of the QPO is best-fit by a small hot region, possibly the inner regions of the accre- 
tion disk, where the ballistic accretion stream impacts onto the disk. The lack of any 
excess power around the QPO frequency in the X-ray power spectrum, created using 
simultaneous lightcurves from RXTE, implies that the QPO is not simply reprocessed 
X-ray variability. 

Key words: accretion, accretion discs - binaries: close - stars: individual: Hercules 
X-l - X-rays: stars 



1 INTRODUCTION 

Hercules X-l/HZ Herculis belongs to the sub-class of 
LMXBs known as X-ray pulsar binaries. Hercules X-l is a 
1.4 Mq neutron star with a 1.24 s spin period. It is in a 
1.7 day binary orbit with it's ~ 2Mq companion star, HZ 
Herculis. Due to the high inclination of the system, i ~ 85° , 
it shows deep eclipses on the orbital period. 

The variability of Hercules X-l has been very well stud- 
ied at all wavelengths. In the X-rays the spin of the neutron 
star is observed, modulated by the orbital eclipses and a 
third 35-day cycle. This 35-day X-ray periodicity was first 
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discovered by Tananbaum et ah (1972) using the UHURU 
satelli te as on X-ray on-state followed by an extended off- 
state. Jones & Forman (1976) later reported a weaker short- 
on state during the off-state. The 35-day periodicity consists 
of a main-on state, lasting about 10 days, and a short-on 
state lasting 5 days, separa ted by periods where the X-ray 
flux drops to almost zero (Scott & Leahy 1999). The X- 
ray pulsations are clearly seen in the X-rays during the on- 
states but appear with a much lower amplitude during the 
off states, as expected if the cycle is caused by the obscura- 
tion of the central regions of the accretion disk. The evolu- 
tion of the pulse profile between and during these on-states 
has been well studied but the origins and evolution of it's 
complex structure remain poorly defined. 

While the existence and form of the 35-day cycle is well 
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documented, the mechanism responsible is less well under- 
stood. It is thought to be caused by a warped, precessing 
accretio n disk that occu lts the neutron star at certain disk 
phases. Pringle (1996) showed that a strongly irradiated 
accretion disk is unstable to warping. The radiation pres- 
sure acting on the inner regions of the accretion disk forms 
a warped precessing accretio n disk as seen in Hercu les X- 
l(|Maloney & Begelman 1997j; [Wijers fc Pringle 1999| ). One 



other model is that the warp i n the disk can be formed by 
an X-ray driven coronal wind ( Schandl fc Meyer 1994 ). In 



this model, the X-rays from the central compact object irra- 
diate the inner regions of the accretion disk and produce a 
hot corona, when the sound velocity of the gas exceeds the 
escape velocity of the potential well the material leaves the 
binar y as a coronal wind, exert ing a repulsive force on the 
disk. Schandl & Meyer (1994) showed that, with suitable 



boundary conditions, they could create a suitable warp to 
explain the general features of the 35-day cycle. This model 
can also explain the observed anomalous and pre-eclipse dips 
as an interaction of the accretion stream with the warped 
disk, causing disk thickness increases that temporarily ob- 
scure the line o f sight to the X-ray source, causing a dip in 
the X-ray flux flgchandl 199rj ). 

At optical wavelengths, much of the light is due to X-ray 
heating of the inner face of the companion star and other 
regions of the binary. This effect is seen to modulate on 
the binary phase and with the 35-day super-period, imply- 
ing that the companion star is irradiated even during the 
observed X-ray off states. The most notable evidence for re- 
processing of X-rays on the companion star of Hercules X-l 
comes from the observed optical pulsations (Middleditch & 
Nelson |l976j ). 

We have found a series of low frequency optical Quasi- 
Periodic Oscillation (QPO) with a central frequency of 35 



mHz. These QP Os have also been observed in the UV (Boro- 
son et il. 2000) with a frequency of 45 mHz. While these 
observations are not simultaneous with the optical data pre- 
sented here, the similarity of the frequencies implies a similar 
production mechanism. QPOs have been observed in many 
binary systems, most noticeably in the X-ray observations of 
both persistent and transient X-ray binaries. In Z-sources, 
such as Scorpius X-l and Cygnus X-2, the QPOs are seen 
to change in both frequency and power as the mass transfer 
rate onto the compact object changes ( Hasinger fc van dei 
Klis 19^|). In transient sources similar X-ray QPOs are seen. 
The transient X-ray pulsar, XTE J1858+034 sh ows 110 mHz 
QPOs with an RMS amplitude of 3.7-7.8 %. |Chakrabarty 
(1998)|]found both X-ray and optical QPOs in the X-ray 
pulsar 4U 1626-67, with a frequency of 48 mHz. The opti- 
cal QPO in this system is most likely due to reprocessing 
off the innerface of the strongly irradiated companion star 
of an X-ray QPO. This X-ray QPO is possibly caused by 
'blobs' on Kepleri an orbits outside the corotation radius of 
the neutron star (Kommers, Chakrabarty & Lewin 1998). 

In this paper we describe the optical characteristics of 
the QPO, in order to constrain the position of the region 
in the binary producing the variability and infer possible 
production mechanisms. 



Parameter 


Start 


End 


UT 


07:40 


08:50 


MJD 


51000.32018 


51000.36879 


Binary phase 


0.407 


0.435 


35-day phase 


0.524 


0.526 



Table 1. Summary of optical obse rvations. The binary phase is 



calculated using the ephemeris of Deet er et al. (1991 



35— day phase uses the ephemeris from ( 3till et al. 2001 ) 



and the 



2 OBSERVATIONS 

The optical data was taken on 1998 July 6, using the Low 
Resolution Imaging Spectrograph (LRIS; Oke et al. (1995)) 



on the 10-m Keck II telescope on Mauna Kea, Hawaii. A 
summary of the observations is given in Table 1 The LRIS 
was used with a 5.2 arcsecond slit masked with aluminized 
mylar tape to form a square aperture. The 300/5000 grating 
used has a mean dispersion of 2.55 A/pixel in the range 
3600A - 9200A. We used a novel data acquisition system to 
obtain more than 58,000 spectra of Hercules X-l, in the form 
of a continuous byte stream. The mean integration time was 
measured to be 72.075ms and there is no dead-time between 
individual spectra. In addition to the 2048-pixel spectrum, 
25-pixel under-scan and 75-pixel over-scan regions were used 
to measure the CCD bias level, which we subtracted from 
each spectrum. The noise for a given pixel was calculated 
using a readout noise of 6.3 e~ and a gain of 4.7 e~/ADU. 
Cosmic rays were rejected with a threshold of 10-sigma from 
the de-biased frames. A master flatfield image was created 
by finding the median of 700 individual flatfield spectra. 
This image showed no deviations above 0.3 % in all but 3 
pixels. It was decided that it was therefore not necessary 
to flatfield individual spectra. Calibration arc spectra and 
spectra of the sky in the region of the object were taken at 
regular intervals. 

The background sky spectrum, which accounts for ~ 
0.9% of the total flux was subtracted using the best fit sky 
spectrum. Sky spectra were taken at the beginning and end 
of the run, so that long timescale variations could be de- 
tected. The mean and variable components of the sky spec- 
trum were found by creating a lightcurve for each pixel and 
extracting the mean and gradient of this lightcurve. These 
coefficients were then filtered in wavelength, with a running 
median filter of width 101 pixels. The gradient of the data 
was found to be almost zero, as the sky brightness didn't 
change noticeably during the run. The arc calibration was 
done by fitting a second order polynomial to 7 lines in a me- 
dian spectra of exposures of Hg and Ar lamps. Arc spectra 
from the beginning and end of the run were used to take into 
account any drifts in the wavelength scale. The wavelength 
calibration was applied using the MOLLY spectral analysis 
package. The individual spectra were flux cal ibrated us ing 
exposures taken of the standard star, Feige 56 (Oke 1990). A 
low order polynomial fit was found to the median of all the 
individual flux star spectra and this calibration was again 
applied using MOLLY. 
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2.1 Time calibration 

Due to the nature of the continuous stream of data, it 
was impossible to attach accurate timemarks to each in- 
dividual spectrum. Individual time marks, accurate to ap- 
proximately 200-ms, were placed after every other spectrum 
using the computer clock. In order to establish an abso- 
lute time reference, secondary timestamps were obtained on 
many occasions during the 5-night observing run. To create 
these timestamps, an incandescent lamp that illuminated 
the CCD was modulated in synch with radio broadcast time 
signals. From these we found an ephemeris and mean expo- 
sure time. These were checked against the computer's time 
marks and found to be in excellent agreement. 

In order to increase the signal to noise of the individual 
spectra and make the data-set more manageable, 7 individ- 
ual spectra were binned together to give 9484 binned spectra 
with a time resolution of 0.5045 seconds. 



2.2 X-ray observations 

The observations took place on 1998 July 6, using the 
PCA onboard RXTE. Lightcurves were created from the 
E_250us_128M_0_lS mode data-set. The energy range of 
the PCA is 2-100 keV. The lightcurves were corrected to the 
solar-system barycentre, using the subroutine fxbary in the 
FTOOLS software suite. They were also corrected for the mo- 
tion of the neutron star in its binary orbit around its 2.2 M^ 
companion H Z Herculis, using the orbital ephemeris of 3til] 
et al. (2001) 



3 RESULTS 

The continuum lightcurve created from our optical dataset, 
in the wavelength range 5000A - 5800A, is shown in the 
upper panel of Figure 1. A subset of this lightcurve, clearly 
showing the QPO is shown in the lower panel of Figure 1. 

We used two methods, one temporal and one spectral, to 
characterise the variability of the source. We created power 
spectra of the source to determine the temporal character- 
istics of the variability, then we used this information to de- 
termine the spectral characteristics by means of a variability 
spectrum. 



3.1 The power spectrum 

We created power spectra for the total data-set by simple 
sine-curve fitting to the input lightcurve, in a method anal- 



ogou s to the Lomb-Scargle periodogram (Press & Rybick: 



198S). However in our method, optimal weights were applied 



when fitting each sine-curve. This method was preferred to 
fast Fourier transforms (FFTs) due to the presence of irregu- 
lar gaps in the data. The resulting power spectrum is shown 
in Figure ^| The power spectrum was represented as a sum 
of three terms; (1) white noise power at all frequencies, (2) a 
'red noise' power- law describing the intrinsic flickering of the 
source, and (3) a hump-like peaked noise feature, caused by 
the summation of the individual transient QPOs, described 
as a Gaussian. These features can be described numerically 
as, 




UTC 




UTC 

Figure 1. Top panel, lightcurve of the optical variability in Her- 
cules X-l for the entire data-set. Bottom panel, lightcurve of a 
subset of the data clearly showing the QPO. The units of flux 
density are mjy in the wavelength range 5000 - 5800 A. Each 
data point is the average of approximately 7 individual spectra. 
No error bars are shown in the top panel for clarity. 
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Parameter 



unit 



value 



Pwn 

a. 

Ppl 

Pqpo 

fo 
A/ 

Integrated QPO 
Power 



10 -3 mjy 2 Hz~ 



~ 3 mjy 2 Hz~ 
" 3 mjy 2 Hz~ 

mHz 

mHz 
10~ 2 mjy 2 



5.09 
2.25 
12.0 
49.30 
35.2 
17.4 
21.5 



0.14 
0.15 
+0.13 
-0.18 
+ 1.4 
1.4 
10.4 
7.1 
+4.3 
1 

+3.7 
2.6 
+6.4 
4.5 



Table 2. Best-fit values to model for power spectrum of total 
dataset, shown in Equation [l]. The best fit \ 2 was 7410.15 for 
4687 degrees of freedom. 



where P wn is the power of the white noise level, P p i is the 
power at f p , fixed to be 0.1 Hz, of the 'red noise' power-law 
with spectral index a. P qpo is the peak power of the Gaus- 
sian QPO feature with a central frequency fo and standard 
deviation A/. (Note: while a lorentzian is normally used to 
characterize a QPO, we have reverted to using a Gaussian 
to characterize the randomness of the position of the central 
frequencies of the individual QPOs) 
The coherence if the QPO is low, 



Q 



A 

A/ 



(2) 



suggesting that the sinusoidal components change period, 
phase or amplitude substantially in just a few cycles, as can 
be seen in the bottom panel of Figure 1. 

The best fit values of the free parameters are shown in 
Table [| These were found using an amoeba code to search 



through the parameter space, see ^ress et al. (1992) 



for a 

discussion of amoeba codes. The best fit has a x 2 °f 7410.15 
for 4687 degrees of freedom. We determined the 1-parameter 
1-sigma confidence regions for each parameter by fixing a 
parameter and then optimising the fit until the badness of 
fit is equal to Xaim- Xaim i s determined using the value of 



Xmin, where, 



2 2 , \rnin 

A.aim \min > 



N 



(3) 



where N — P is the number of degrees of freedom in the fit. 



3.2 Dynamic power spectrum 

We used the Gabor transform (Heil & Walnut 1989,1990) to 
create dynamic power spectra to investigate the time evo- 
lution of the power spectrum determined from our optical 
data. This is especially valuable due to the low coherency of 
the oscillations, to determine if the low coherency is caused 
by the evolution of the central frequency of the oscilla- 
tion, perhaps caused by material rapidly moving through 
the inner regions of the accretion disk, or whether the low 
coherency is caused by randomly occurring events based 
around a central frequency, as would be the case for short- 
lived blobs of material on Keplerian orbits within the disk. 

The discrete Gabor transform is a wavelet transform 
that allows us to determine the power spectrum over a sub- 
s et of the data, u sing a Gaussian window on the lightcurve 
( Boyd et al. 1995 ) . We have chosen the 3-sigma width of the 
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Figure 4. The X-ray power spectrum for our RXTE observations 
simultaneous with the optical observations. The power spectrum 
has been normalised, following the method of Leahy et al. (1983) 
and 



the Poisson level has been removed. The power spectrum 
shows no evidence for excess power in the region of the optical 
QPO. However an excess of power is clearly seen at 0.01 Hz. 



window to be 750-s. The data set contained 9484 data points, 
which have been filtered with a running median filter with 
a filter width of 101-s to remove the slowly varying compo- 
nent. For computational purposes, the resulting lightcurve 
was padded with 6900 randomly distributed points with the 
same mean and variance as the data set. The results of this 
analysis are shown in Figure ^, where the main image shows 
the results of the Gabor transform, the lower panel shows 
the integrated power at a given time and the right-hand 
panel shows the average power spectrum. 

The dynamic power spectrum shows transient QPOs in 
the frequency range of the peaked noise feature. While there 
is little evidence of frequency evolution of the individual 
features, it is possible that the feature that appears at ~ 15 
mHz at T = 2800-s, evolves into the feature at ~ 46 mHz 
at T = 4500-s. From the coherency of the individual QPOS, 
we see that the average lifetime of the features is 200-300-s, 
or ten times the Kepler period of the features. 



3.3 A reprocessed X-ray QPO ? 

In order to determine if the optical QPO is a reprocessed ver- 
sion of an X-ray QPO, we calculated power spectra of the 
X-ray data simultaneous with our optical data. We chopped 
the total X-ray lightcurve into 512-s sectio ns and created in - 
dividual Leahy-normalised power spectra (Leahy et al. 1983) 



for each lightcurve. We then combined the separate power 
spectra to create a total power spectrum for the observa- 
tions. This method is more successful at determining the 
power in the QPO, which spans a broad range of frequen- 
cies and is only coherent on short timescales. 

The average power spectrum, shown in Figure ^, shows 
no evidence for the QPO at 35 mHz, implying the optical 
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best fit 




1,05 0.1 0.2 

frequency (Hz) 



residuals of fit 




frequency (Hz) 

Figure 2. Top, the binned power spectrum for the total dataset, showing the best fit model (solid line), with a reduced x 2 = 1-581 
for 4687 degrees of freedom. The white noise, power law and QPO components (dashed lines) are also shown based on their best fit 
parameters. Bottom, the residuals of the fit, showing the full resolution of the power spectrum. 



QPO is not simply an X-ray QPO that is reprocessed into an 
optical flux by some region in the binary. However there is 
evidence for some excess power at 0.01 Hz, which is possibly 
the X-ray counterpart of our slow variability. 



3.4 The spectrum of the variability 

We created variability spectrum-lightcurve pairs that, when 
re-combined, reproduce the entire spectral time series, using 



a method similar to that used by Eracleous & Home (1996) 
to decompose the flaring and quiescent spectra of AE Aqr 
The model of the spectral time series, M(A, t), is written as 



M (A, t) = So (A) + Si (A) Li (t) + S 2 (A) L 2 (t) , 



(4) 



where So (A) is the average spectrum, Si (A) and Li (t), are 
the spectrum and lightcurve of the first variable component 
and S 2 (A) and L2 (t) are the spectrum and lightcurve of the 
second variable component. Li (t) and L2 (t) are normalised 
to a mean value of zero and unit variance. This normalisation 
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Figure 3. Results from performing a Gabor transform on the filtered optical dataset. The main panel shows the dynamic power spectrum. 
The right-hand panel shows the average power spectrum and the lower panel shows the integrated power as a function of time. 



ensures that the spectral components Si (A) and 52 (A) have 
units of mjy and do not contribute to the time averaged 
spectrum. 

The model is compared to the observed spectral time 
series and the badness of fit calculated using a \ 2 analysis. 
We solved iteratively for each component of the model, until 
a suitable convergence criterion was met, namely that the 
percentage change in the global % 2 was l ess than 0.1% over 
3 successive iterations. In order to check that a global min- 
imum had been reached, we chose several different sets of 
initial values, each of which converged onto the same final 
values. 



In order to separate the variability into slow and fast 
components we filtered the lig htcurves, using Savitzky- 
Golay filters (Press et al. 1992) of different widths. We 



low-pass filtered the lightcurve L\(t) to remove variations 
faster than 80-s. This component then measures primarily 
the red noise power-law component identified in the power 
spectrum. Similarly we low-pass filtered the lightcurve L2 (i) 
to remove variations faster than 5-s, which are primarily the 
white noise component. This means that Li (t) responds to 
variations between 80-s and 5-s, which include the tail of the 
red noise and the QPO. 

The results of the the fit to the spectral time series are 
shown in Figure ||. In the first variability lightcurve, L\ (t) 
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the fluctuations seen with frequencies around 10 mHz are 
resolved. The spectrum of this variability, Si (A), is bluer 
than the average spectrum, So (A). The rms amplitude of 
this variable component increases from 1% of the mean flux 
at the red end of the spectrum to 2% at the blue end. The 
second variability lightcurve, L2 (t), shows the variability of 
the 35 mHz QPO. The features are transient, showing the 
QPO is not coherent over more than a few cycles. The spec- 
trum of this variable component, S2 (A), has an rms ampli- 
tude of 0.5 % at the red end of the spectrum to 1 % in the 
blue. The spectral features in the two variability spectra are 
very similar. The Balmer lines appear to be narrower than 
in the average spectrum, implying they have variable wings, 
while the cores of the lines vary less than the average spec- 
trum. The Nil + CH blend at 4641 A also appears to show 
a strong variable component. This line is interesting as it is 



a blen d of two lines formed b y Bowen fluorescence (Schac- 



ter, Fil jpenko & Kahn 198£ ) and is thought to be caused 
by irradiation of the binary by the central X-ray source. In 
contrast the He II line at 4685 A is not apparent in either 
variability spectrum, implying that it is not highly variable. 

These results agree well with the power spectra analy- 
sis in the previous section. The average flux in the second 
variability lightcurve, shown in the third panel in Figure |B| 
is in good agreement with the flux in the QPO and the vari- 
ability lightcurve shows variability on the same timescales 
as the QPO period. We conclude from this that the second 
variability spectrum, S2 (A), represents the spectrum of the 
QPO. 



3.5 Spectral fits to variability spectra 

The next stage of the analysis was to fit model spectra to 
the data to fit physical parameters to the observed average 
and variability spectra from the previous analysis. To do this 
we used two models to see which best represents the data. 
Both models use blackbody spectra to describe the observed 
continuum fluxes. 

(i) A single blackbody: 

This will tell us the temperature and size of the emitting 
region. The model spectrum is defined at a wavelength A as, 



F„(A) 



ttR 
~D 



2 B V (X,T) 



(5) 



for a region with a radius, R, a temperature, T, at a distance 
D. This model is most useful for a region that has a roughly 
constant temperature for the duration of the observations. 

(ii) The difference of two blackbodies: 

The second model is fitting the difference of two black- 
bodies to the data. This involves fitting the difference of two 
blackbodies with different temperatures, Ti and T2, where 
T2 > Ti. The model spectrum is defined at a wavelength A 
as, 



F„(A) 



nR 2 



[B„(A,T 2 ) -B„(A,Ti)] 



(6) 



for a region with a radius, R, at a distance D. This model 
is most useful for describing a region whose temperature 
varies during the observations. These variations are limited 
to coming from the same region, ie a single region where the 
temperature varies between Ti and T2. 



These models were fitted to regions of the spectra that 
are least contaminated with lines. The best fit values, to- 
gether with their 1-parameter 1-sigma confidence regions are 
shown in Table ^| The data together with the best fit mod- 
els and the regions used are shown in Figure ^| The value 
of x 2 is large in all three cases, due to the high signal-to- 
noise of each spectrum revealing intrinsic deviations from a 
Blackbody. However, with the confidence intervals re-scaled 
as shown in Equation ^, the parameters remain reasonably 
well constrained. 



4 DISCUSSION 

We have studied the X-ray and optical variability in the ac- 
creting pulsar Hercules X-l and found evidence for a series 
of optical QPOs that comprise a broad peaked noise feature 
in the overall power spectrum. While the coherence of this 
feature is low, Q ~ 2, the individual QPOs that comprise the 
peak have much higher coherencies, as seen in the lightcurve 
in the lower panel of Figure 1 and the widths of the features 
in the dynamic power spectrum. We have characterised this 
noise feature as having a Gaussian shape in the power spec- 
trum. This has been fitted along with the white noise level 
and a power law to describe the variability of Hercules X-l 
at low frequencies. The best fit to this model has a peak 
frequency of 35.2 mHz. 

We have extracted the spectrum of the variability from 
the frequency range of the QPO. We have also determined 
the average spectrum and the spectrum of any slower vari- 
ability, that has characteristic frequencies from the red noise, 
power-law region of the power spectrum, as shown in Fig- 
ure ^. The variability spectra are bluer and fainter than the 
average spectrum, indicating they are from small, high tem- 
perature regions. The average spectrum is predominantly 
that of the companion star, HZ Herculis. It shows strong 
absorption in the Balmer lines, as well as some He I absorp- 
tion lines. There is some evidence of the irradiation of the 
companion in the form of a weak He II (A4686) emission line 
and Bowen fluorescence lines (A4634, A4641, A4642), which 
can be used as a diagnostic of the physical conditions in the 
binary (Bchacter, Filppenko & Kahn 1989). The variability 
spectra both show a bluer continuum than the average spec- 
trum with Balmer and Helium I absorption lines, however 
in these spectra the wings of the lines show more variability 
than the cores. The Bowen lines also show more variability 
than the average spectrum, as expected if the variability is 
coming from a region that has been heated by the irradiation 
of the X-ray source. 

Our simple blackbody fits to the observed spectra, while 
not particularly physical, enable us to constrain models for 
their production. The best fit model to the average spectrum 
is a single temperature blackbody, with a mean temperature 
of 16500 K. The size of the emitting region is 8.6 x 10 22 cm 2 . 
While this temperature is inconsistent with the temperature 
of the companion star in the absence of irradiation, which 
has been observed during eclipse of the accretion disk to 
have a temperature of ~ 7000 K, it is most likely to be due 
to the strongly irradiated inner face of the companion star. 
The size of emitting region would cover approximately 70% 
of the innerface of the companion star, assuming a distance 
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Figure 5. Top panels, the average spectrum, So (A), (left) and its associated lightcurve (right), which is set to be unity. Centre panels, 
the first variability spectrum, Si (A), (left) and its associated lightcurve, L\ (t), (right). Bottom panels, the same as above for the second 
component of the variability. Typical error bars are shown on the spectra. The first variability lightcurve is detrended and filtered with 
a width 40-s. The second variability lightcurve is filtered with a width of 5-s. 



of 5 kpc and an effective Roche lobe radius of 3.5 x 10 11 cm 



(Eggleton 198^) 



The best fit model to the variability spectra of the 
slowly varying component also enable us to place constraints 
on the size and possible locations of the emitting region. The 
slowly varying component is best fit with the derivative of 
a blackbody model, although the single temperature black- 
body is almost as good a fit. In this model the source varies 
between 7600 K and 31000 K, from a region of 3.6 x 10 20 
cm 2 . One candidate for this component is the same region 
on the companion star as the average spectrum, but with 



a smaller area, possibly the most highly irradiated regions 
close to the inner Lagrangian point. 

The best fit model to the spectrum of the QPO is that 
of a single temperature blackbody with a mean temperature 
of ~ 33600 K, with an area of 1.6 x 10 20 cm 2 . The data is 
not consistent with the model of a derivative of a blackbody 
as used to describe the slowly varying component. The fit 
converges to almost the same solution as with the single 
temperature Blackbody. This leads us to consider whether 
or not the QPO is produced in a physically distinct region 
from that of the slowly varying component. The QPO pro- 
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unit 


average spectrum 
So 


slow variability 
Si 


fast variability 

s 2 


model 1: 
Temp 
Radius 

x 2 /ioo 


K 

10 10 cm 


16500 ±i°S 
29.36 ±°:^ 
9016 


37000 + 1000 

3/UUU _ 1000 

1 64 + - 04 
L -"* -0.03 

3.919 


33600 t\T 
1 26 + 03 

J -" ) -0.04 

2.966 


model 2: 
Ti 
T 2 
Radius 

X 2 /99 


K 
K 

10 10 cm 


600 +19 °° 

16500 tZ 

29 4 + - 2 
za '^ -0.2 

9107 


7600 + 9900 

31400 + 5700 
cti^uu _ 14600 

1 89 + mf 
i.oa _ 24 

3.904 


600 + 44 ?° 

— mf 

33600 

1 26 + ' 04 
± - ZD -0.04 

2.996 



Table 3. Summary of best fit parameters to the spectral fitting of the spectra from the variability model, together with their 1-parameter 
1-sigma confidence regions. A distance of 5 kpc was assumed throughout. 



during region is a small, high temperature region, probably 
in the accretion disk surrounding the compact object. There- 
fore if the frequency of the QPO is the Kepler frequency of 
material orbiting around the neutron star, the Kepler radius 
of the material orbiting with a frequency of 35 mHz around 
a 1.3 Mq neutron star would be 3.3 x 10 9 cm. This would 
appear to be too small to contain the region responsible for 
the observed spectrum. However if we take the range of fre- 
quencies included in the QPO, we find that the Kepler radii 
covered are 9.4 x 10 s - 5 x 10 10 cm, which is large enough 
to encompass the observed region. 

The existence of a range of preferred azimuths in the 
disk explains the observations of the QPO frequencies. There 
are two mechanisms whereby this anisotropy could cause 
the observed QPO. The first is the material in the preferred 
region of the disk is orbiting the neutron star and we are 
seeing the irradiated innerface of the region, as it passes 
across the backside of the disk. The second explanation is 
that the azimuths are fixed in the plane of the binary and 
that inhomogeneities in the material traveling through these 
regions gives rise to enhanced emission. This material would 
orbit the neutron star and emit when it passes through the 
preferred azimuths. 

Schandl (1996) modelled the warp in the accretion disk 
to explain the occurance of the anomalous and pre-eclipse 
dips, using a coronal wind model. In this model it is the vari- 
ation with synodic period of the penetration depth of the 
stream into the inner regions of the accretion that accounts 
for an increase in material along the line of sight to the bi- 
nary that accounts for the observed dips. The Keck II obser- 
vations were taken at the synodic orbital phase of 0.894. At 
this phase the accretion stream can penetrate down to the 
inner regions of the accretion disk. According to Schandl's 
calculations at this synodic period the stream will penetrate 
the disk at a radius of ~ 5.1 x 10 10 cm. This is in excellent 
agreement with the lowest Kepler frequency of the observed 
QPO. This model predicts that the thickness of the disk will 
increase as the stream interacts with the disk. If this energy 
is dissipated on a thermal timescale, then for a thin disk, 
tdyn ~ ottth, so that a ~ Q~ ■ Thus our observation that 
the coherence of an individual feature in the dynamic power 
spectrum last for ~ 10 x tdyn naturally leads to a value of 
a ~ 0.1. 

Such a model, where the disk-stream interaction point 



leads to a blob of material orbiting in the disk, close to 
the corotation radius, could also explain the observations 
of the QPO in 4U 1626-67. In this system, the QPO has a 
similarly l ow coherence, but is observable over m ore than 
a decade (Kommers, Chakrabarty & Lewin 1998), which is 



easily explained if the accretion stream causes a relatively 
stable preferred radius in the disk. However in Hercules X-l 
we do not observe the QPO in the X-ray lightcurves, as in 
4U 1626-67, although this could be due to the viewing angle 
during the X-ray off state. Observations during the X-ray on 
state are necessary, with simultaneous optical observations 
to determine if the optical QPO is indeed a reprocessed X- 
ray QPO. 
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